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Abstract 

Starting  with  Picht's  solution  of  the  wave  equation  in  the  form 
of  an  integral  of  plane  waves  over  a  caustic  we  derive  expressions  for 
the  solution  wMch  depend  entirely  on  the  geometry  of  the  problem.     An 
asymptotic  evaluation  of  Picht's  integral  for   high  frequencies   gives 
the   desired  result  which  can  have  a   number   of  different  forms  depending 
on  the  precise  region  where  the  field  is  being  observed.      In  particular, 
a  change  in  the  asymptotic  field  occurs  in  going  from  the  regular  ray 
field  to  the  caustic   or  from  a  regular  part  of  the  caustic   to  a  cusp. 
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1.     Introduction 

In  this  paper  we  shall  consider  a  particular  class  of 
solutions  of  the  wave  equation  in  two  dimensions, 

2     2 

(1,1)     1^  +  i_^  +  kVu  =0  (n  =  constant), 

in  a  form  which  is  especially  suited  for  dealing  with  the  case  where 
the  wavelength  is  asymptotically  small.  It  is  well  known  that  in 
the  small  wavelength  limit  of  diffraction  or  propagation  phenomena 
the  geometrical  features  dominate  the  underlying  physical  situation. 
Accordingly,  the  solutions  we  shall  use  are  completely  determined  by 
geometrical  considerations. 

The  geometry  of  a  field  can  be  described  optically  by  a  set  of 
rays  which  may  be  regarded  as  the  paths  of  energy  flow.  VJhere  the 
rays  converge  to  a  focus,  or  form  an  en'/elope  (also  known  as  a  caustic), 
it  may  be  assumed  that  the  energy  density  or  the  field  amplitude  is 
large.  Where  there  is  a  region  containing  no  rays  at  all,  the  energy 
does  not  propagate.  Given  the  ray  system,  the  field  amplitude  in  the 
small  wavelength  limit  can  be  computed  at  most  points  of  space  by 
well  known  methods.  In  fact,  the  optical  field  amplitude  has  been 
calculated  even  at  points  on  the  smooth  part  of  the  caustic   , 
although  the  field  becomes  infinite  in  the  small  wavelength  limit  on 
a  caustic.  It  is  known  that  the  dependence  of  the  asymptotic  field 
on  the  wavelength  changes  suddenly  along  a  ray  as  it  passes  through 
a  point  on  a  caustic,  and  a  change  must  occur  again  in  a  different  way 
where  the  caustic  has  a  cusp.  We  shall' investigate  the  behavior  of 
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tlie  asymptotic  field  in  those  sin,^lar  regions  by  means  of  a  class 

[2] 

of  explicitly  f^iven  fields  whose   construction  is  due  to  J.   Picht  ^  •■ . 

r3i 

The  Picht  fields  are  a  generalization  of  one  used  by  BebyB  ^  ••  for 
investigating  the  case  of  a  ray  system  which  converges  to  a  point 
focus. 

If  any  field  vrhatever  has  a  "geometrical  optics  limit  where  the 
ray  system  is  regular,   there  i-iill  exist  a   superposition  of  Picht  fields 
which,   together,    have  the  same   asymptotic  limit.     Moreover,  the  Picht 
field  will  have  the  same  asyinptotic  limit  en  the  caustics  of  the  ray 
system  as  the  given  field,   a  fact  vrhich  can  be  easily  seen  from  an 
application  of  Green's  theorem*-  ■' .     Thus,    the  conclusions  we  shall 
reach  are  actually  more  general  th^n  might  be  expected  at  first. 

What  we  desire  here  is  a  purely  geometrical  description  of  the 
field,   a  description  which  in  principle  would  allow  us  to  find  the 
field  amplitudes  graphically  from  a  graph.ic  reproduction  of  the 
corrsspondiny  ray  system, 

2.     The  Picht  Solution 

To  construct  the  exact  solution  of  (l.l)   given  by  Picht,  wo  must 
first  determine  the  ray  system  of  the  corresponding  geometrical optics 
problem.     This  can  be  done  whether  vre  are  interested  in  diffraction 
by  an  obstacle  in  a  homogeneous  medium,  in  which  case  the  rays  are 
straight  lines,   or  in  propagation  through  an  inhomogcneous  medium,   in 
which  case  the  rays  are  curved  lines.     In  any  event  the  ray  system  will 
generally  form  one  or  more  envelopes  known  as  caustics.     The  caustic 
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curves,  on  the  other  hand,  completely  determine  the  ray  system,  for 
at  each  point  on  a  caustic  there  is  a  tangent  ray,  and  every  ray  is 
tangent  to  one  of  the  caustics. 

The  geometrical  optics  apprordmation  fails  at  the  csustics 
because  there,  the  true  wave  amplitude  grows  indefinitely  vjith  the 
frequency.  Our  attention  will  be  fi:ced  upon  the  immediate  neighborhood 
of  tlio  caustic.  Thus,  if  the  rays  of  the  system  are  curved,  we 
can  replace  them  "by   the  ray  system  composed  of  straight  line  rays 
having  the  same  envelopes  as  the  original  system;  for  inf initesim^al 
distances  there  is  no  distinction  between  straight  lines  and  sufficiently 
smooth  curved  linos. 

Let  P  be  the  position  vector  to   the  point  at  which  the  field  is 
to  be  observedj  let  R(s)  be  the  radius  vector  to  the  caustic  C,  and 
s  the  arclength  parameter  along  C.  Let  p(s)  be  the  radius  of  curvature 
of  C,  and  lot  A(s)  be  a  quantity  determined  by  the  initial  distribution 
of  the  field  amplitude,  e.g.  at  the  scurce.  V7e  have,  for  the  unit 
tangent  vector  to  C, 

(2.1)     ^=T(3)        (fTl^^l). 

The  solution  of  Picht,  corresponding  to  the  caustic  C,  is  then 


(2.2)   u(P)  =    A(s). 


exp  ik  [(?-R(s)).f(s)+s]'  ^^ 


pTsT  , 


where  we  have  taken  the  index  of  refraction  n  to  be  one  for  convenience. 


(For  Figure  1,   see  page  19a.) 


The  function  exp  ik  P  .  T  ,  where  T  is  an  arbitrarj'  constant 
VGctor  of  unit  amplitude  is  a  plane  wave  solution  of  the  wave  equation 
(1.1),     It  follows  that    the  integrand  in  (2.2)  and  thus  tliat  u(?) 
itself  is  a  solution  of  (1.3).     Because  of  the  special  form  of  the 
integrand  in  (2.2)  u(?)  asyrrptota  cally  approeches  the  appropriate 
geometrical  optics  field  at   the  point  P  as  k  becomes  large. 

3.     The  A53^ptotic  Solution 

The  Picht  solution  (2.2)  is  valid  everywhere  for  large  v.--lues 
of  k  insofar  as   the  caustic  C  determines  tlie  field  u(rj.     In  what 
follows  we  shall  investigate  in  more  detail  the  asymptotic 
evaluation  of   (2.2)   for  large  k  so  that  W3  may  obtain  approximate 
expressions  for  u(P),  which  are  simpler  and  depend  more  explicitly 
upon  the  ray  geometiy  near  P^ than  upon  (2.2).     In  the  discussion  we 
shall  distinguish  between  three  kinds  of  regions  in  which  P  can  be 
located.     The  point  P  is   in  the    "  lit"  region  if  a  real  ray  tangent 
to  C  at  a  real  point  passes   through  P.      If  no  such  ray  fix^m  C  can 
reach  P,  ire  say  that  P  is   in  the    "  shadow"      region.     Finally,    P  can 
lie  on  the  caustic  C  itself. 

In  the  lit  region  we  can  evaluate  the  integral  in  (2.2)  by 
stationary'  phase  without  running  into  complications.     The  derivative 
with  respect  to  s  of  the  phase  of  the   integrand  in  (2.2)   is 


(3.1)  [P^-  R(3)]   .  If(s)/p(s) 


where  N(s)  is  the  unit  normal  to  the  caustic  C.  Clearly  (3.1)  will 

vanish  for  s  =  s  if  the  vector  P  -  R(s^)  is  parallel  to  the  unit 
o  o 
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vector  "^Cs  )   tangent  to  C,   i.e.,   when  P  -  H(s^)  lies  along  the  ray- 
through  P.     Since  P  is  in  the  lit  region  this  is  possible  for  real 
values  of  the  vectors,   and  the   corresponding  real  value  s^  of  the  arc 
length  parameter  s  is  the  stationary  point.     Let  us  assume  that   F 
is  on  a  ray  betxreen  the  source  and  the  stationary  point  s     on  the 

The  second  derivative  of  the  phase  or  the  first  derivative  of 

(3.1)  at  the  stationary  point  s     is 

(3.2)  -D(s^)/p^(s^)       , 

where  D(s   )  is  the  Tncgnitude  of  the  vector  P  -  R(s   )• 

By  expanding  the  phase  in  a  power  series  about  the  stationary 
point  s  ,   and  keeping  only-  teims  up  to  the  quadratic,  we  obtain  the 
approximation. 

(3.3)  u(F)  roj     A(s)  expJik  rD(s^)  +  s^-D(s^)(s-s^)V2p^(s^)    |ds/p(s). 

The  important  contribution  from  (30)  comes  from  the  immediate 
neighborhood  of  s  .  Thus,  we  can  take  the  interval  of  integration 
from  -00  to  +00  in  s  without  much  effect  on  the  result.  Now  we  define 


3. 


9-9^=1  ds/p(s)/vj(s -s^)/p(s^) 


so  that  ds/p(s)  =  d&,  whence 
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(3.h) 


u(P)  <'^|gxp  ik   p(s^)+s^ 


fOO 


J 


M&)  expl  -ikD(s^)(&-e^)^/2U9 


■  cc 


'•^A(9   )  \  exp  Ik  p(3   )  +s 


Ik  p(3^) 


rOO 


exp  J  - 


ikD(s  )e'^/2 

0 


.  de. 


-00 


Finally,   we  obtain  frorr.  (B.'i)    (see  Appendix  l): 

(3.5)  u(?)   rOy2n/kD(s^)  A(s^)  exp  J  ik|D(s^)   +  s^     -  inA  I    . 

The  quantity  D(3   )    can  be  recognized  as  the  radius  of  curvature 
of  the  wavefront  through  P  at  the   ray  which  is   tangent  to  the   caustic 
at  3   J   for  the  caustic  is  alvjays  the  locus  of  the  centers  of  curvature 
of  any  wavefront.     The  equation  of  the  iravefix^nt  through  P  is 


(3.6)  D(3)  +  s  =  D(s^)  +  s^     : 


here  D(s)  is  the  distance  clong  the  ray  at  s  from  the  caustic  to  the 
wavefax^nt,   i.e.,   the  radius  of  curvature  of  the  wavefront  at  the 
wavefront  end  of  the  ray  through  s. 

If  we  set 


(3.6)  A(3^)   =  A(9^)  /ZJT^x  exp(iTiA)         , 


then  we  have 


(3.7) 


u('P)   =  A(9^)  D~^'^^(Sq)   expJik b(s^)  +5^  I 


The  normalized  amplitude  A(&   )  will  be  determined  by  the  angular 
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distribution  of  energy  at  the  source.  We  recognize  (3.7)  as  the 
geometrical  optics  solution  for  the  lit  region;  in  particular,  note 
that  the  anplitude  varies  inversely  as  the  square  root  of  the  radius 
of  curvature  of  the  wavofront. 

If  we  assume  that  P  is  in  the  lit  region,  not  between  the  source 
and  tte  caustic,  but  rather  in  such  a  position  that  the  stationary 
point  on  the  caustic  lies  on  a  ray  between  the  source  and  P,  the 
results  at  each  stage  of  the  calculation  are  the  same  except  for  the 
sign  of  tlie  quadrr^tic  term  of  the  power  series  of  the  phase  in  (3.3). 
The  effect  is  accounted  for  by  multiplying  (3.5)  or  (3.7)  by 
i  =  exp(in/2).  Thus,  we  have  the  well  known  jump  in  phase  of  the 
geometrxcal  optics  solution  when  a  ray  passes  through  the  caustic. 

VJhen  the  point  P  is  in  the  shadow  region  of  the  caustic,  there 
is  no  real  value  of  s  on  the  caustic  C  for  which  the  vector  P  -  R(.s; 
is  parallel  to  the  tangent  T(s)j  i.e.,  one  cannot  construct  a  real 
tangent  to  C  which  passes  through  P.  However,  if  the  radius  vector 
R(s)  of  C  is  an  analytic  function  of  s,  it  can  be  continued  into  the 
complex  plane  analytically.  Then,  defining 

(3.8)    ^1   =  t(s)/p(s) 

by  analytic  continuation  for  complex  values  of  s,  in  general  one  can  find 
a  complex  value  s  of  s  such  that  the  phase  is  stationary,  i.e., 


|V-l(s„)]. 


(3.9)    |F-"?(s„)l.N(Sj,)  -0 
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The   contour  C  must  be  deformed  until  it  passes  through  s   .     The 
power  series  approximation  for  the  phase  up  to  quadratic  terms  can 
again  be  used,   and  the  resulting  integral  evaluated  as  before.     The 
approximation  will  resemble   (3.7)   except  that,   owing  to  the   complex 
phase,  an  exponential  decay  factor  will  be  present.     Thus,   in  general 
the  solution  decays  exponents  a] ly  in  the  shadow  region. 

The  quantities  involved  in  this   case  are  geometrical  only  in  an 
extended  sense  by  analytic  continuation.  If  the  caustic  C  is  an 
arbitrarj'  curve  and  P  is  in  the  neighborhood  of  C,   in  the  sh'dow 
region,   then  we  may  obtain  an  approximate  evaluation  of  (2.2)   in 
terms  of  real  geometrical  quantities.      This  can  be  done  by  replacing 
C  by  its  circle  of  curvature  at  the  point  on  C  nearest  the  point  P 
and  then  evaluating  (2.2)  as  if  the  integration  contour  were  the 
circle  of  curveture.     The  evaluation  of  (2.2)   for  the  case  of  a 
circular  caustic  will  be  carried  out  later. 

When  P  is  at  a  regular  point  s     of  the  caustic  C,    s     is  also  a 

0  '     o 

stationary  point,    for  obviously  (3.1)  vanishes  in  this  case.     However, 
since  D     is  zero,    (3«2)  must  vanish  also.     To  obtain  the  asjTnptotic 
evaluation  of  (2.2)  we  must,   therefore,  use  the  next  non-vanishing 
term  of  the  power  series  of  the  phase  in  (2.2).     The  third  derivative 
of  the  phase  is 


(3 


.no)  -r^(g/(p*l/p3^   (?-f).N.3(f/p3)(P-R).?*l/p^  . 


At  s     we  have 

0 


P  -  R  =  0   ) 


-  2D  - 


thus,    {3»2D)  becomss  l/p     and 


r°°  /-  '\ 


•  ds/p(s)    . 


(3.11)     u(?)^A(s   )(exp  iks^)  e:zp  \    ik(s  -  s^)^/6p^(s^) 

O  0       J 

-00 

Proceeding  as  before,  we  obtain  from  (3.11) 

;oo  V 

exp  Jik  p(s   )  e^/6l    d9 
-ex.        ^  ^ 

or 

^         2p  (^)cos   n/6  ,  ,  ,       * 

u{?)oj 3 ---*y  a(s   )   3r.p  ik  s  (See  Appendix  1)     * 

3[kp{sJ/SY^        "^ 

If  the  caustic  C  is  "  flat"  (i.e.,  p  =  oo)  at  the  point  s  ,  so 
that  the  third  dorivative  (3.10)  vanishes,  we  must  use  the  next  order 
term  in  the  power  series  expansion  of  the  phase  in  (2.2).  For  the 
general  case  in  which  the  first  n  terms  of  the  power  series  vanish 
we  use  the  (n  +  l)   term,  and  the  procedure  is  the  same  as  tliat  for 
the  more  regular  cases  just  described. 

Where  the  caustic  C  has  a  cusp  the  radius  vector  R(s)  has  a 
singularity.  It  is  convenient  to  assume  that  the  origin  of  our 
coordinate  system  is  at  the  cusp  which  we  can  take  to  be  the  point  s  =  0 
on  C.  Then  the  cusp  is  at  R(C )  =0.  We  can  also  assume  that  the 
coordinate  system  is  so  oriented  that  the  tangent  vector  T(s)  has  a 
horizontal  limiting  position  at  s  =0j  i.e.,  at  s  =0 

T(0)  =  (1,0)   . 

Geometrically  a  cusp  is  a  point  at  which  two  branches  of  the 
curve  C  meet.  The  field  should  be  computed  separately  for  each  branch 
and  the  total  field  obtained  by  superposing  the  two  contributions. 


This  result  agrees  with  that  on  page  160  of  [l]  . 
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Thus,  we  consider  now  the  contribution  of  a  single  branch  for 
which  s  ranges  from  zero  through  all  positive  values. 

It  is  appropriate  to  assime  a  singularity  in  the  ci^rvature 
l/p(s)  at  the  cusp  md  we  shall  ass'jme  that  the  singularity  is  algebraic! 

(3.13)    l/p(s)  =  \s  +  tis*^  +  a  function  of  higher  order  in  s. 

Here  p  >  a  and  we  also  assume  that  a  >  -1,  which  is  necessary  in 
order  that  T(s)  exist  at  s  =  0. 

Now  we  can  write  in  general  for  the  components  of  R(s): 

(3. DO    t{s)   =  J  I   cos  J  dt/p(t)  da",  sin  J  dt/p(t)  d&  \   , 

'o      ^  '{j     '^  J 

whence 

(3.1?)    7(s)  -I    cos  J  dt/p(t)  ,  sin  I  dt/p(t)  (   , 

and 

(3.16)    ^   =  N(3)/p(3)  =^  -sin  I  dt/p(t)  ,  cos  j  dt/p(t)^  /p(s) 


From  (3.IU)  and  (3.15),  using  the  power  series  expansions  for  the  sin 
and  cos  functions  about  zero  we  obtain 

(3.17)  f  -  Is  -  p/2(a+l)^(2a-.3)V°''^  -  ...  , 

r\/(a-H)(a+2)'j  s*"*^  *  rn/(p+l)(p+2)ls^^^  +  ...I 

(3.18)  f  =  h  -ir|x/(a+l)j  s''*^  Jn/(P-Hl)j  sH^  +  ...  , 

[x/(a+l)1  s'"'^  +  L/(P+1)1  s^*^  ^  ...1  . 
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Writing  P  =  (x,y)  we  then  have  for  the  phase 

(3.19)        [P -r].T  +s  =  X  +   fyX/(a+l)]s''''^  +  [yvL/(P+l)]s^*^ 

.[xX^x/(c.l)(p.l)]   s^-^P-^2^[xV(a^2)(2a.3)]   s^""'^ 

\Ihen  P  is  at  the  cusp  x  =  y  =  0,   and  the  phase  becomes 

2a+3 


+ 


(3.20) 


rxV(a+2)(a+3)] 


Along  a  ray  y  =  0,  and  the  phase  becomes  near  the  cusp 


(3.21)    X  -  lxXV2(a+l)^  s^""^  +  ... 


-  rxxV2(a-M)2J 


We  can  also  consider  various  other  cases  in  which  we  approach  the  cusp 
along  different  paths,  e.g.,  along  x  =  0,  and  in  each  case  we  take 
from  (3.19)  only  the  lowest  order  term  in  3.  For  these  other  approaches 
a  more  precise  relationship  between  a  and  p  must  be  given  to  decide 
which  term  is  actually  of  lowest  order. 

In  (2.2)  we  can  replace  A(s)  by  A(0)  and  the  phase  by  (3.20), 

(3.21)  or  by  a  special  case  of  (3.19),  depending  on  the  particular 
limiting  case  vie   wish  to  compute.  Consider,  for  exam.ple,  the  case  of 
X  =  y  =  0.  From  (2.2),  (3.20)  and  Appendix  1  we  have  : 

(3.22)  u(P)^A(C)(sgn  X)\x\^^^'^^^\2a.3y^-'^^^^^''*^^ 

.j^-(a.l)/(2a.3}(^,2)^"'^^/^2"-^3)pr(^^l)/(2a.3)] 
.  exp(in(a+l)/2(2a+3)|  . 
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The  contribution  of  the  other  branch  of  the   cusp  in  the  c?3e  of 
coinplete  symmetry  about  the  horizontal  is  calculated  in  the  same  way. 
However,  wo  must  be  careful  that  the  integration  is   in  such  a 
direction  that  T(s)  turns  so  that   its  anj^le  with  the  horizontal 


9  = 


ds/p(s) 


increases.  On  the  other  hand,  if  we  assume  that  the  second  branch  is 
also  given  for  s,  going  from  zero  through  positive  values,  it  is 
clear  geometrically  that  p(s)  must  have  a  sign  opposite  that  of  p(s) 
on  the  first  branch.  The  integration  over  C  for  the  second  branch 
then  must  be  taken  in  the  direction,  -ds.  Now  to  first  order  in  s  we 
have  from  (3.13) 

sgn  X  =  sgn  p(s)  . 

Thus,  in  the  case  x  =  y  =  0  an  examination  of  (3.22)  indicates  that  a 
sign  change  occurs  for  integration  in  the  direction  ds,  and,  therefore, 
integration  in  the  direction  -ds  cancels  the  effect  of  the  sign  change 
exactly.  The  contribution  of  the  second  branch  then  will  be  exactly'' 
the  same  as  that  of  the  first  branch,  and  the  total  field  at  the  cusp 
will  be  exactly  twice  the  value  given  by  (3.22)  with  sgn  X  equal  to  +. 

It  is  interesting  to  consider  also  the  case  where  P  is  on  the 
ray  y  =  0  through  the  cusp.  From  (2.2),  (3.21)  and  Appendix  1. 


7 


(3.23)       u(F)  rNJ  (sgn  X)A(0)(Tt/2k|x|  )^  exp. 


ikx  -  (iii/li)sgn  x 
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The  value  of  u(P)  given  by  the  expression  (3.23)  when  sgn  X  is  + 
is  exactly  one-h^lf  that  given  by  (3.5)  for  a  point  P  on  an  ordinary 
ray  in  the  lit  region.     This  factor  of  one-lialf  occurs  because  we 
are  considering  just  the  effect  of  one  branch  of  the  cusp.     The  other 
branch  of  the  cusp  produces  a  contribution  to  u(P)  identical  with  that 
given  by  (3.23).     The   sum  of  the  two   contributions  will  then  result  in 
a  field  value  for  u(P)  which  is  the  same  as  that  occurring  for  an 
ordinarj""  position  of  P  in  a  lit  region. 

For  the  case  of  P  on  the  line,  x  =  0  orthogonal  to  the  ray  y  =  0, 
we  choose  the  second  and  third  terms  on  the  right  side  of  (3.19)  for 
the  Dhases 

(3.2li)       [yX/(a+l)]  s'^-'^  +  LaP+d]  s^^^       . 


From  (2.2)  and  (3.2li)  we  have 

■  •  ds 


(3.2?)       u(P)/N^XA(0)      I     s*  exp  ik   lyXAa+l)  3°*^+    W(P+1)  s^*^ 


k  ry\Aa-^l)]s"*^*  [y^Ap+1)]^ 


=  fu(0)Aa+l)]  f  expiJcl[y\Aa*l)]t'+LAM)]'t^^*^^^''''^|dtf. 

For  certain  relations  between  a  and  p,  e.g.,  for  (p+l)Aa+l)  an  integer, 
(3.25)  can  be  evaluated  easily  by  stationary  phase.  For  the 
contribution  of  the  other  branch  of  the  cusp  we  multiply  the  result 
(3,25)  by  -1  and  replace  X  and  |j.  by  -X  and  -\i. 
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II.     An  Example;  A   Circular  Caustic  ard  a  Point  Focus 

Let  us  consider  a  simple  example  of  a  caustic  vrlth  no   cusp. 
Let  C  bo  a  circle  of  radius  p  vrith  its   center  at   the  origin,     Th3n  C 
is  given  by 

(U.l)         R(s)   =f  p  CGsCs/p),   p  sin  (s/p)  j  , 
Also  we  liave: 

T  =  j -sin(s/p),    cos(3/p)j 

N  =  -lcos(s/p),   sin(s/p)j    , 

and  p  is  the  constant  radius  of  curvature  of  E.     The  relation  (2.2) 

becomes 

u(P)   =    j     A(s)   exp^ikl-x  sinCs/p) +  y  cos(s/p)+  si  y  ds/p 

(U.2)  ^^ 


3(»)   exp 


where 


<ik  -X  3in(3/p)  +  y  cos(s/p)+3    > 

q)jik     r  sin  i<^  -9)  +pg    I  d&     , 

A  s(9)     =  3(©),   r  =   |f|   and  (/  is  the  direction  angle  of  F, 
and  s  =  p&.     Tlie  angle  9     for  which  P  -  R  is  normal  to  IJ  can  be 
computed  quite  easily  or  seen  at  once  geometrically: 

(h.3)  »o  "  ^  ■  cos"^(p/r)   . 

The  distance  D     is  easily  seen  to  be 

°  1 

(h.]4)  D^  =  (r^-P^)2     . 

.^ 

From  (3*5)  for  a  point  P  in  the  lit  region  r  >  p  i-re  have 

(U.5)     u(?)-^J2TtA(r2-p2)2-    B(&^)explik  (r%2)^  +  p((i(-cos"^(p/r))-i:iA 


where  we  have  used  p&  =  s     and  ()i.3). 
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The  result   (h.S)  holds  only  in  the  case  that  a  single 
portion  of  the  wavefront  corresponriing  to  C  passes   through  P. 
In  general,   a   spiral  shaped  wavefront  winds  around  C  and,   as  it 
propagates,   one  turn  of  the   spiral  after  another  will  pass  through 
P.     At  each  passage  the  stationary  point  is   such  that  the 
corresponding  angle  9     is  increased  by  2tt.     This  fact  is  predictable 
by  the  observation  that  the  term  cos'7p/r)of  (h.3)   is  ambiguous. 
Thus,    in  general,   the  field   at  P  is 


li        .  1 


(U.6)     n(?)rjy    B(9^+2Tm)l2nA(r%^)^     expJik  j(r^-p^)'^+p(0-cos~-^/ty2nn)J| 

n=-oo 

Of  course  BfS-  +2Tin),  which  is  given  by  the  source,  may  vanish  for  an 
o 

argument  sufficiently  large  in  magnitude  and  reduce  (Ii.6)  to  a  finite 
sum. 


For  P  on  C  where  r  =  p  we  have  from  (3.12) 

1. 
(U.7)    u(P)^|; 


2p(l/3)(cos   n/6)/3[kp/6pl  A(p{!f)exp(iic(:0) 


for  a  single  turn  of  the  spiral  vravefront.     In  general, 

1       00  _ 

(U.8)     u("p)^|2p(l/3)(co3   TT/6)/3r!:p/6pj^  A(p(^+2nTip)Gxp  ik(piZf+2Tin)    . 

In  the  case  of  a  point  focus  the  radius  p  of  C  becomes  zero, 
and  we  might  consider  this  as  a  limiting  case  of  a  circular  caustic. 
Setting  p  =  0  in  (li.2)  we  obtain 


(U.9)  u(P)r^j  B(9)  exp  Ilk  r  sin(0-&)  I  d©  . 

C 
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Now  in  the  limit  r  ->  0  a  stationary  point  6  will  lag  behind  0 

by  n/2.  It  is  convenient  to  replace  the  function  B(&)  in  (h.?) 

by  an  angular  source  distribution  runction  which  determines  the 

proportion  of  energy  assigned  to  each  ray.  Since  the  direction  of 

a  ray  through  P  is  p,   this  czn  be  accomplished  by  a  change  of  variable 

(h.lO)         9  +  n/2  =  C  . 

Correspondingly  we  set 

(U.ll)      F(4)  =  BU-n/2)   . 

Then  (U.9)  becomes 

(h.l2)  u(P)  =   I     F(©)exp  Jikr  cos(0-e)f     d©     , 


u(P)  =   I     F(©)explikr  cos(0-e)| 

-71  V 


vrhere  we  liave  assumed  only  a  single  integration  about  C  which  implies 
that  F(6)   is  periodic  cf  j^eriod  2vi.      If  this  vrere  not  true  we  should 
be  forced  to  extend  the   integral  (U.3  2)  over  a  larger  interval.     We 
can  check  on  the  meaning  of  F(&)  with  respect  to  our  source  by 

evaluating  (U.12)  directly  by  stationary  phase: 

1 
-^  -  (  \ 

(U.13)  u(P)  r^(2nAr)''  7(v^)exp  •  ^t  +  InM     , 

v:hich  holds  when  P  is  not  near  the  focus. 

The  result  (b.12)  agrees  with  the  Debye  solution  of  t'ne  perfect 

r3i 

focus  problem'^-'. 

In  case  p  is  not  zero  and  the  point  r  is  inside  the  circular 
caustic  we  have 

(U.lh)         r  <  p  , 
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and  the  stationary  point  given  by  (h.3)  is  complex.     Since 


cos(9^-9   }   =  -2 
or 


it  is  easily  found  that  in  the  complex  case 


(U.i?) 


&   =  ,d  +  i  l0£ 
o  ^  ' 


(p-Vp^  )  /r|  . 


Then  (I4.5)  becomes  in  this  case: 


(U.16)         u(P)  ro 


2nA(p'^-r")2jB(6^) 

1  _  -N 

.  exp|-k|(p^-r^r  +  p  log(p/r-V^p-r  /r)     +  iW/ 


where  B(9  )  is  the  analytic  continuation  of  the  source  distribution 
to  the  complex  angle  8  given  by  (U.l5). 


5.  Concliisions 

It  can  be  seen  from  an  examination  of  the  various  fomulas  for 
the  field  amplitudes  given  in  Sections  3  and  h   that  the  field  is 
completely  specified  in  the  small  vfavelength  limit  by  the  local 
geometry  determined  by  the  caustics  of  the  ray  system.  For  points 
in  the  lit  region  or  on  a  smooth  part  of  a  caustic  it  is  necessary  to 
know  the  local  geometry  only  to  the  ord.er  of  the  radius  of  curvature 
of  the  caustic.  In  the  neighborhood  of  a  cusp  of  a  caustic  the 
necessary  information  must  include  the  order  of  contact  of  the  two 
caustic  branches.  In  the  shadow  region  an  approximate  value  for  the 
field  amplitude  near  the  caustic  is  again  given  in  terms  of  the 
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curvature  of  the   caustic.     For  points  f-jrther  away  in  the  shadow 
region  the  small  wavelength  limit  of  the  field  amplitude  is 
exponentially  small,   and  it  can  be  argued  that  non-zero  values 
for   this  case  occur  only  ulien  the  wave  length  is  sufficiently  long 
to  introduce  non-geometrical  considerations   into  the  problem.     In 
case  a  point  in  the  shadow  region  is  near  enough  to  several  points  of 
closest  approach  on  a   caustic,    the  field  should  consist  of  a 
super-position  of  fields  determined  by  (h.l6)   for  each  such  point 
of  closest  approach, 

(For  Figxires  2  and  3,   see  page  19a.) 

Appendix  I 
For  the  sake  of  completeness  we  consider  the  evaluation  of  an 
integral  of  the  form 

fCD 

(A.l)  I  =         exp  jikas    I  ds  , 

where  a  is  a  real  constant  and  X  is  a  r^ial  positive  constant. 

After  a  change  of  variable 
1     1 
(A.2)  s  =   (|a|k)"^t?^  exp  |(Ti/2X)3gn  aj 

in  (A.l)  we  obtain 

1       ,  ..    /-ioo     1 


I  = 
(A.3) 


=  |(|aIk)'^/J|cxpr(n/2\)sgn  ajj  J  f^      exp(-r)  ^X 

=   l(|a|k)"'vJ|expr(Ti/2X)sgn  all    P  (^)      • 
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Figiire  1. 
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Figiire  2. 


Figure  3< 
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